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Abstract The study of spatially explieit integro-difference systems when the 
local population dynamics are given in terms of discrete-time generations mod- 
els, has gained considerable attention over the past two decades. These nonlin- 
ear systems arise naturally in the study of the spatial dispersal of organisms. 
The brunt of the mathematical research on these systems, particularly, when 
dealing with cooperative systems, has focused on the study of the existence 
of traveling wave solutions and the characterization of their spreading speed. 
Here, we characterize the minimum propagation (spreading) speed, via the 
convergence of initial data to wave solutions, for a large class of non cooper- 
ative nonlinear systems of integro-difference equations. The spreading speed 
turns out to be the one with slowest speed from a family of non-constant trav- 
eling wave solutions. The applicability of these theoretical results is illustrated 
through the explicit study of an integro-difference system with local popula- 
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tion dynamics governed by HasscU and Comins' non-cooperative competition 
model (1976). The corresponding integro-difference nonlinear systems that re- 
sults from the redistribution of individuals via a dispersal kernel is shown to 
satisfy conditions that guarantee the existence of minimum speeds and travel- 
ing waves. This article is dedicated to Simon A. Levin whose contributions to 
the fields of ecology, evolutionary biology, and the environmental sciences have 
driven and inspired the research of generations of mathematicians, mathemat- 
ical biologists, and life and social scientists, around the world for over four 
decades. 

Keywords minimum speed • traveling wave • dispersal • biological invasion ■ 
integro-difference systems 

Mathematics Subject Classification (2000) MSG 39 All • 92D25 • 
92D40 



1 Introduction 

Finding and developing macroscopic descriptions for the dynamics and behav- 
ior of heterogeneous large ensembles of individuals subject to ecological forces 
like dispersal continues to provide challenges and opportunities for mathe- 
matical and biological scientists. Over the past century, particular attention 
has been placed on the study of the role played by dispersal in shaping plant 
communities, in helping understand biological invasions, in assisting in the 
quantification and control of the spread of infectious disease, or in disentan- 
gling the dynamics of marine open-ocean and intertidal systems, to name but 
a few examples. The work of pioneers like Aronson [1,2], Fisher [7], Hadcler 
[8-10], Kolmogorov [15], Levin [18], Okubo [19], Skellam [29], Slobodkin [14], 
Weinberger [33] and the subsequent cadre of distinguished mathematicians 
and theoreticians across the world who have worked at this interface, set not 
only the foundation of an important and fertile area of interdisciplinary re- 
search (ecology, mathematics, and evolutionary biology) but in the process 
it has inspired novel mathematical research while being re-energized by un- 
solved questions in emerging fields like urban ecology and sustainability and 
the challenges and opportunities posed by the growing body of research on the 
co-evolving dynamics of socio-biological systems [5] [17]. 

The theme of finding mathematical macroscopic descriptions for the spatial 
dynamics of heterogeneous large-ensembles of populations was set in "motion" 
by the fundamental ecological contributions of Skellam (1951) [29], Kierstad 
and Slobodkin (1953) [14], Levin and Paine (1974) [18], Okubo (1980) [19], 
and others. The study of integro-difference equations dispersal models in the 
mathematical literature has its origins in the study of the coupled spatial 
dynamics of organisms with discrete primarily non-overlapping (but see [30]) 
local dynamics with dispersal processes modeled via re-distribution kernels [1, 
2]. 
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Early models for the dispersal of invasive species used nonlinear reaction- 
diffusion equations, with the prototype provided by Fisher's Equation [7]. The 
primary motivation or emphasis have always been on the characterization of 
the speed of propagation of species invading unoccupied habitats. The seminal 
contributions of Fisher [7] and Kolmogorov, Petrowski, and Piscounov [15], 
Aronson and Weinberger [1,2] jointly handled the mathematical challenges 
posed by their efforts to classify rigoroiisly the concept of speed of propaga- 
tion for quite general continuous in space and discrete time, integro-difference 
equations. Weinberger [33] and Lui [24] research expanded the mathematical 
foundation for the theory of spreading speeds and traveling waves, through 
their analysis of traveling waves via the convergence of initial data to wave 
solutions, in the context of cooperative operators. Recently, Weinberger, Lewis 
and Li made additional contributions [31,20,21,32]. The mathematical analy- 
ses of integro-difference spatially explicit systems enhances the understanding 
of the dynamics of introduced species like weeds or pests in terrestrial sys- 
tems, or the study of the impact of dominant alien species in freshwatcrs 
while generating additional challenges and opportunities to mathematicians, 
whose interests, are driven by the study of challenging dynamical systems. 

The pervasiveness of overcompensation in biological systems implies that 
integro-difference equations models are in general non-cooperative, and there- 
fore, existing theoretical work has yet to address effectively the mathematical 
consequences of non-cooperative local dynamics on dispersal. In other words, 
the incorporation of biological forces /mechanisms that drive population over- 
compensation leads to mathematical models whose dynamics have yet to be 
satisfactorily teased out in the context of relevant biological settings. Deep 
mathematical challenges remain [16]. The research in this manuscript does 
not start in the vacuum since relevant mathematical work for non-cooperative 
systems has been carried out by several researchers. Thieme [27] showed, in 
the context of a general model with non-monotone growth functions, that the 
asymptotic spreading speed could still be obtained with the aid of carefully 
constructed monotone functions. Hsu and Zhao [13] and Li, Lewis and Wein- 
berger [22] just extended the theory of spreading speeds in the context of non- 
monotone integro-difference equations. Their extensions relied on two meth- 
ods: the construction of two monotone operators (with appropriate properties) 
and the application of fixed point theorems in Banach spaces-an approach also 
used in Ma [26] and Wang [28] to establish the existence of traveling wave so- 
lutions of reaction-diffusion equations. The results in this manuscript on the 
speed of propagation for non-cooperative systems in the context of integro- 
difference equations rely on the spreading results for monotone systems in 
Weinberger et al. [31]. 

We highlight our results in the context of a two-dimensional nonlinear dis- 
crete system describing the local nonlinear dynamics of two competing species 
with discrete reproduction cycles [11]. The model focuses on the growth and 
spread of these competing species with their population densities at genera- 
tion n and spatial location x being tracked by the state variables Xn{x) and 
Yn{x), respectively. The system is a natural extension of the classical single 
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population "scramble" competition model of Ricker [3] . Specifically, the non- 
spatial interference-competition model of Hassell and Comins is given by the 
following system of coupled nonlinear difference equations: 

where ri, r2, cri, (72 are all positive constants. 

The possibility that individuals in the above two populations may dis- 
perse to different sites is modeled with a redistribution kernel ki{y). Hence, 
a discrete-time model, where individuals interact locally according to Model 
(1), can be naturally formulated via a system of coupled nonlinear integro- 
difference equations. Hence, we have that 

JR 

where the dispersal of the i-species is modeled by a redistribution kernel ki,i = 
1, 2 that depends just on the signed distance x — y, connecting the "birth" y 
location and the "settlement" location x. In other words, ki{y) is a homogenous 
"probability" kernel that satisfies ki{y)dy = 1. 

Since the above system is non-cooperative in general, it is in such a context 
that new results will be formulated, and illustrated but first, we introduce 
the notation that will be used for the explicit mathematical formulation of 
the dynamics of two-interacting, dispersing, and competing populations. Con- 
sequently, (3, jj"^ ,F,F^ ,r,u,v are used to denote vectors in or iV- vector 
valued functions while x,y,^ are used to denote variables in R. The use of 
u — (m*) and v ~ (w*) G allow us to define u > v whenever > 
for all i; and u :$> v whenever > for all i. We further define for any 
r = (r') > 0,r e R^ the i?^-interval 

[0,r] = {u:0<u< r,u € R^} C R^ 

and 

Cr = {u = (u\ u^) : u' G C(R, R), < u'{x) < r\x gR, i = 1, N}, 

where C(R, R) is the set of all continuous functions from R to R. Our focus 
will be on the set C^+, /3+ > 0. 

Specifically, we consider the system of integro-difference equations 

= Q[F{un)]; (3) 

where u„ = «) G C^+, F{u) = (/i(w)); 

Q[F{u)] = {Q'[F{u)]y, 



(1) 



(2) 
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Q'[F{u)]{x)= I h{x-y)fi{u{y))dy; 

Un{x) is the density of individuals at point x and time/generation n; F{u) is 
the density-dependent fecundity (local growth rate) ; and h {x — y) (dispersal 
kernel) models the dispersal of u, assumed to depend only on the signed dis- 
tance x—y between the location of "birth" y and the "settlement" or "landing" 
location x. As noted before, ki{x — y) can be viewed as a probability kernel 
since ki{x)dx = 1. The notation Q[F{un)\ is slightly different from these 
used in [13, 22, 24, 31] and hence those wishing to compare results must account 
for this. Specifically, no F can be found in the standard literature notation. 
Here, F has been included to carry out the proofs involving non- monotone 
systems effectively. 

The integro-difference system (3) models the reproduction and dispersal of 
a time-synchronized species where all individuals first undergo reproduction, 
then redistribute their offspring, and then proceed to reproduce again. The 
goal is to carry out the characterization of the spreading speed in a system 
involving a rather general non-cooperative system (3) as the slowest speed of 
a family of non-constant traveling wave solutions of (3). 



2 Non-cooperative Systems' Results 

The focus is on the characterization of the speeds of propagation for (3) when 
the system is non-cooperative. As it is typical in mathematics, we make use of 
prior results established results for cooperative systems ([31]). The existence 
of two additional monotone operators F^ with the properties that the first 
lies above F and the second below F is required by our method of proof. 
The use of this approach is motivated by the work on non-monotone equa- 
tions carried out in [27,13,22,26,35,28]. Specifically, we observe that F^ can 
be "constructed" via piecewise functions made up of "pieces' of F and the 
incorporation of appropriate constants. If F happens to be monotone, then 
F^ = F. We introduce additional technical assumptions below. The assump- 
tions are critical because the feasibility of the mathematical analysis depends 
on whether or not the components of our problem meet them: 

(HI) For i = 1, ...,N, kiir) > is integrable on R, kiir) = ki{-T),T G M, and 

/jj ki{T)dT = 1, /jj ki{T)e^'^dT < +oo, for all A > 0. 
(H2) (i) Given that F : [0, /J"*"] [0, /S"*"] is a continuous, twice piecewise con- 
tinuous differentiable function with 

« /3- = {(5-) <p= (A) <P+ = (/5+), 

it is assumed that there exist continuous, twice piecewise continuous 
differentiable functions F± = (/f ) : [0,/3+] ^ [0,/3+] such that for 
[0,/3+], 

F-{u) < F{u) < F+{u). 
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(ii) F{0) — 0, F{13) = (3 and there is no other positive equilibrium of Q[F] 
between and /3 (that is, there is no constant v ^ (3 such that F{v) = 
v,Q '^v < P). F^{<S) = 0, F^{P^) = and there is no other positive 
cquihbrium of Q[-F^] between and P'^. F has a finite number of 
equilibria in [0,/3+]. 

(iii) F"^ are nondecreasing functions on [0, and F"^ [u) and F{u) have 
the same Jacobian at 0. 

Assumptions (H1-H2) do not suffice if the goal is to characterize the speeds 
of propagation for (3). The assumption (H3), which includes the requirement 
that the operator grows less than its linearization along the particular func- 
tion Ufj_e~^^, is essential and implies that the operator Q does not display an 
Allee effect for this particular function (see [31]). Assumption (H3), explicitly 
formulated below, is satisfied by several biological systems of interest this will 
be highlighted in our example. Assumption (H3), therefore does not severely 
handicap the usefulness of the results in this manuscript. 

The need for Frobenius' theorem stating that any nonzero irreducible ma- 
trix with nonnegative entries has a unique principal positive eigenvalue with 
a corresponding principal eigenvector "made up" of strictly positive coordi- 
nates is implicit in Assumption (H3) . The formulation of (H3) depends on the 
concept of irreducibility. A matrix is irreducible if it is not similar to a lower 
triangular matrix with two blocks via a permutation (See [12,31]). By reorder- 
ing the coordinates, one can put any matrix into a block lower triangular form, 
then we say that the matrix is in Frobenius form if all the diagonal blocks are 
irreducible (an irreducible matrix consists of the single diagonal block which is 
the matrix itself). Here we use the definition of Frobenius form in Weinberger 
et al. [31]. Following the approach in [31], that for each > 0, the N x N 
matrix that results from the linearization of (3) at 0, namely 

= (bif) = {djfm [ h{s)e^^ds), (4) 
Jr 

where is the {i,j) entry of the matrix, is in Frobenius form ([31]). In other 
words, it is assume that the required reordering has been done for ([31]). 
If we now let A(/i) denote the principal eigenvalue of the first diagonal blocks, 
we reach the formulation required by Assumption (H3): 

(H3) (i) Assume that is in Frobenius form and that the principal eigenvalue, 
A(/x), of the first diagonal block is strictly larger than the principal 
eigenvalues of other diagonal blocks. Further, let's assume that B^ has 
a positive eigenvector i/^ = (t'^) » corresponding to A(/i) with the 
additional requirement that A(0) > 1. 
(ii) For each /x > and a > 0, we let v"^ = (vf) = {min{Pf, uj^a}), and 
assume that 



F^{v^) < Bov^. 
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(iii) For every sufficiently larger positive integer k, there is a small constant 
vector w = (w*) » such that 

F^{u) > {l-^)Bou, uG[0,u], 

It follows from (HI) that A(/u) is an even function. In fact, it was shown by 

Lui [24] that \nX{p) is a convex function and therefore, \nX{p) achieves its 
minimum at /x = and, therefore the assumption that A(0) > 1 implies that 
lnA(/u) > 0. The statement in Proposition 1 below which is critical to the rest 
of analysis that leads to the main result and it involves the following function 
of the largest principal eigenvalue A(/l() 

^{li) = -lnA(/x) >0. 
M 

Part (5) of Proposition 1 highlights the use of this function in the construction 
of lower solutions and estimates of the traveling wave solutions. 

Proposition 1 Assume that (HI) — {H3) hold. Then 

1 ^(/i) — > oo as /i — > 0; 

2 is decreasing as fi = 0+; 

3 changes sign at most once on (0,oo) 

4 ^{p) has a minimum c* > 0. 

5 For each a > c* , there exist Ac > and 7 G (1, 2) such that 

<P(ylc) = c, <P(7/lc) < c. 

Parts (l)-(4) of Proposition 1 are essentially due to Lui [24]. However, Lui's 
results only guarantee that c* > 0. The proof of the strict inequality, that is, 
that c* > 0, is found in the Appendix. Since A(/i) is a simple root of the 
characteristic equation of an irreducible block, it can be shown that A(/u) is 
twice continuously differentiable on K. Part (5) is a direct consequence of the 
results stated in Parts (l)-(4). 

A traveling wave solution m„ of (3) is defined as a solution of the form 
UniO = u{£_—cn),u £ C(R, R^). The theorems that guarantee the existence of 
traveling wave solutions for cooperative systems have been already established 
(e.g. [21]). It also has been established that the asymptotic spreading speed, 
for such systems, can be characterized as the speed of the slowest non-constant 
traveling wave solution for monotone operators [21] and for scalar equations 
[34,33,13,22]. 

We start with the statement of Theorem 1, the main theorem, which gener- 
alizes results previously established for cooperative systems to non-cooperative 
systems. Some parts of Theorem 1 such as the asymptotic behavior of travel- 
ing waves are new even for cooperative systems. The new information about 
cooperative systems are, in fact, required to be able to carry out the proofs of 
the results for non-cooperative systems. The details associated with the proof 
of the main result are collected in a series of lemmas and theorems all collected 
in the following sections. 
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The two major new contributions in this paper are: 1) for a large class of 
non- monotone systems (3), the question of the existence of the minimum speed 
of propagation is settled (Theorem l(i-ii)) and this speed is characterize as the 
speed of the slowest non-constant traveling wave solution (Theorem l(iii-v)); 
and 2) in the case of competition model, a direct application of the main theo- 
rem helps identify simple and meaningful conditions needed for the existence 
of traveling waves with the minimum speed of propagation (Theorem 3). That 
is, what is required to guarantee the success of a biological invasion. It is worth 
re-iterating that the application of the results in this manuscript to the study 
of relevant monotone operators case [21,33] does give additional information. 
In fact, these results help explicitly characterize the asymptotic behavior of 
traveling waves via the careful analysis of eigenvalues and upper-lower solu- 
tions. This analysis was not done before [21,33] most likely because the focus 
was exclusively in establishing the existence of traveling waves. The results and 
analysis for the n- dimensional case is typically harder. Our approach works 
because the analysis of the n dimensional case is closely related to structure of 
the eigenvalues and corresponding eigenvectors, an analysis that is embedded 
in our study of the relevant monotone operators. 

The following theorem summarizes the main results. 

Theorem 1 Assume (HI) — {H3) hold. Then the following statements are 
valid: 

(i) For any uq € with compact support and < uq /?, the solution Un of 
(3) satisfies 

lim sup Un{x) = 0, for c > c* 

\x\>nc 

(a) For any strictly positive vector w € ]S.^ , there is a positive Roj with the 
property that if uq £ and uq>uj on an interval of length 2R^, then the 
solution Un{x) of (3) satisfies 

j3~ < liminf inf Un{x) < for < c < c* 

(Hi) For each c> c* (3) admits a traveling wave solution u{^—cn) = (w*(^— cn)) 
such that < w(0 < /?+, ^ e R, 

(3~ < liminf m(^) < limsupu(^) < (3~^ 
limj^cxa u{^) = and 

lim u(Oe'^=^ =z^Ae. (5) 

//, in addition, F is non-decreasing on Cp, then u is non-increasing on M. 
(iv) For c = c* (3) admits a non-constant traveling wave solution — cn) = 

{u'{^ - cn)) such that < m(C) < C e M. 
(v) For < c < c* (3) does not admit a traveling wave solution u„(5) = 
— cn) such that u € with liminf5_,._oo w(^) » and u{-\-<x>) = 0. 
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Remark 1 When F is monotone, = F,I3^ = p. 

Remark 2 The assumption that F has a finite number of equihbria in [0, ,9+] 
is only used in the proof of Theorem 1 (iv), and can be further relaxed. In 
fact, as long as for some component i and a sufficiently small positive number 
5,u= > with = 5 are not equilibria of F, the conclusion is still valid 
from the proof. 

We shall establish Theorem 1 in Sections 3 and 4. 



3 Spreading Speeds 

Our results on the speed of propagation for non-cooperative systems make use 
of Theorem 2 below which collects the properties of the spreading speed c* for 
monotone systems as established in Weinberger, Lewis and Li [31]. Theorem 
2 extends the related spreading results in Lui [24] to systems of monotone 
recursive operators with more than two equilibria. The operator at the center 
of this manuscript may support more than two equilibria with one lying at the 
boundary as in [31] (see Section 5). 

Theorem 2 (Weinberger, Lewis and Li [31] [Lemma 2.2, Theorem 3.1]) As- 
sume [HI) — {H3) hold. Further assume that p{x),i = 1,...,N is non- 
decreasing. Then the following statements are valid: 

(i) For any hq g with compact support and < uq <C /3, the solution Un{x) 
of (3) satisfies 

lim sup Un{x) = 0, for c> c* 

\x\>nc 

(ii) For any strictly positive vector w £ M.^ , there is a positive R^j with the 
property that if uq G Cj3 and uq > uj on an interval of length 2Ru, then the 
solution Un{x) of (3) satisfies 

liminf inf u„(x) — /3, for < c < c* 

n^oo \x\<nc 

It is clear that Q[-F='=] are monotone (order preserving) on That is, if 

u,v G 0^+ and u{x) < v{x),x € M, then 

Q[F±(u)](a;) < Q[F'^{v)]{x), x eR. 
Further, for u= (u') G and x gR, we have 

f-{u{x))<f,{u{x))<f+{u{x)),i = l,...,N. 

and therefore 

Q[F-{u)]{x) < Q[F{u)]{x) < Q[F+{u)]{x), xgR. 
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We arc now able to establish Part (i) and (ii) of Theorem 1 by following 
essentially the proof for the scalar cases found in [13,22]. 
Proof of Parts (i) and (ii) of Theorem 1. 

Part (i). For a given uq G with compact support, let u„ be the n-th 
iteration of Q[F] starting from uq and let u+ be the n-th iteration of Q[-F~''] 
starting from Uq. By (H2), we have 

< u„{x) < u+(a;),a; e M,n > 0. 

Thus for any c > c* , it follows from Theorem 2 (i) that 

lim sup u^{x) = 0, 

\x\>nc 

and hence 

lim sup |u„(a;)| = 0, 

\x\>nc 

Part (ii). Let Un,u^ be the n-th iteration of Q[-F],Q[-F+] starting from uq 
respectively. Let Vq = minjuQ, i = 1, ...,N. Then vq = (vq) G C^- . Letting 

u~ denote the n-th iteration of starting from vq and observing that 

vo < uo and (3~ </?< /?"'", from (H2), we have that 

u~{x) < Un{x) < u^{x),x e K, n > 0. 

Theorem 2 (ii) states that for any strictly positive constant w, there is a 
positive (choose the larger one between the for F+ and the for 
F^) with the property that if uq > w on an interval of length 2i?„. Hence, it 
follows that the solutions u^{x) satisfy 

lim inf inf (x) = , for < c < c* . 

t->oo \x\<tc 

Thus for any c < c*, it follows from Theorem 2 (ii) that 

lim inf inf u^{x) = /3^, 

n-foo |x|<nc 

and consequently, that 

P~ < lim inf inf Un{x) < . 

n^oo \x\>nc 



□ 
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4 Characterization of c* as the slowest speeds of traveling waves 

A non-constant solution of (3) is a traveling wave of speed c provided that 

it has the form Un{x) = u{x — cn), where u G C(]R,R^) and, of course, if 
it satisfies (3). By substituting this form into (3), it follows that u{^) must 
satisfy the following system of equations. 

uio = Qc[F{u)m = mnuM)) ■■= Qinnm + c) (6) 

In this section we complete the proof of Theorem 1 (iii), (iv) and (v), that 
is, the portion of our main result that characterizes the spread speed c* as 
the speed of the slowest member of a family of non-constant traveling wave 
solutions. This is an extension of prior results for monotone operators [21] and 
for scalar equations [34,33,13,22]. 



4.1 Upper and lower solutions 

In this subsection, we shall verify that f/)"*" and defined below are the upper 
and lower solutions of (6) respectively. These solutions are only continuous on 
M. Upper and lower solutions of this type have been frequently used in the 
literature (see Diekmann [6], Weinberger [34], Lui [24], Weinberger, Lewis and 
Li [31], Rass and RadclifFe [25], Weng and Zhao [36], more recently by Ma [26] 
and Wang [28]). In particular, the explicit use of upper vector- valued solutions 
can be traced to the work in [24,31,25,36]; for lower vector-valued solutions, 
in the context of multi-type epidemic models, to the work in [25]; and in [36] 
in the context of multi-type SIS epidemic models. Our construction of ^+ and 
(t)~, the upper and lower solutions of (6), is motivated by the research in these 
references. 

Our verification of the lower and upper solutions for n- dimensional systems 
is new and different from the above mentioned references. The details follow 
below. 

Let c>c*, l<7<2,g>l and recall the definitions of Ac and 
as utilized in Proposition 1. The corresponding positive eigenvectors ^-nd 
u^Aa of for the eigenvalues when n = Ac, 'yAc can therefore be identified. 

Define 

where 

</.+ =min{A,!^le-^=«}, CgK; 

and 
where 

cpr = max{0, u\^e-^'^ - qv^A^e-^^'^}, ^ G R. 
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In M- In 4^ 

It is clear that if ^ < then (f>+{^) = l3,; and if ^ > then 

<t>tiO - ^^i.e-^^f. Similarly, if ^ < ln(g^)^^3T)X ^^^^ = ^nd 

if i > ln(g^)^^3^ then </.r(0 = iy\^e-^^i - qvl^^^e-^^^^ . 
We choose q> 1 large enough so that 

\n{q^) In 4^ 

> 



(7 - 1)A, -A, 
and therefore 

<^+(O>0r(e),?eK. 

We verify in the two lemmas below that and ^~ are upper and lower 
solutions of (6) rcspcctiv(4>-. It is assumed that Lemma 1 is valid when F is 
monotone. In this case, = F, (3^ = (3. 

Lemma 1 Assume F is monotone and {HI) — {H3) hold. For any c > c* , 
then ^+ is an upper solution of Qc [F] . That is 

In^ 

Proof Let 4* = Then = A if ^ < and </)+(^) = i^^^e-^^? if 

e > Note that <f>+{0 < v\^e-^^i, ^ e M. 
In view of (113) we have, for ^ e R 

TV N 

Thus, for ^ e K, in view of (4), (113), Proposition 1, we obtain that 



(7) 



N 

Q'[F{^'^M + c) < e-^=(«+^) J2 "aA 

= e-Mi+c)e^c^{Ac)^i^^ 
= v^e-^'^^+'^'^e^'^'' 

— ^Ac^ 

On the other hand, since (t>'l{£,) < I3i,i = 1, N , we have for ^ e M 

2lF(0+)](e + c)<A. (8) 

Thus, we have for ^ e M 

Qim^M) = Q'm^m +c)< ct>t{o- (q) 
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This completes the proof of Lemma 1. 

□ 

In order to verify the lower solution, the following estimate for F is needed. 
For N = 1,2, Lemma 2 can be found in [28]. 

Lemma 2 Assume {HI — Hi) hold. There exist positive constants Di,i = 
1, N such that 

N N 

fi{u) > J2 djfi{0)u^ - A u = (ui), u e [0, i = 1, N. 

Proof In a sufficiently small neighborhood of the origin, since F is twice 
continuously differentiable. From the Taylor's Theorem for multi- variable func- 
tions (the big Oh notation version), for u sufficiently small. 

N N 

Mu) =J29jM0W +0{J2{un'), u = {u^),ue[0,f3+],i = l,...,N. 
There exist small e > and -D ■ > such that for < ^ 

N N 

Mu) > J2 9jM0W - D[ " = u € [0, /?+],» = 1, TV. 

For u e [0,/3] and > e, noting that fi{u),J2j'=i9jfi{0)u^ are 

bounded, we always choose a sufficiently large constant D" > such that 

Thus if we let Z), = max{D^, -D-'}, then Lemma 2 is proved. □ 

Lemma 3 Assume (HI) — {H3) hold. For any c > c* if q (which is indepen- 
dent of ^) and that it is sufficiently large, is a lower solution of Qc[F]. 
That is 

Proof Again let = In(g^) (^_\)^^ ■ Hence if ^ < 4* then = 0; 

while if C > 4* ^^^^ 'PiiO = '^Ae^"^"^ ~ Q'^^A^^''^^"^- ^^^y see that 

!^ie-^=« > r (0 > '^le-^^f - g4^^e-^^=€, ^eR,i = 1,...,N. (10) 
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For ^ G R, in view of Lemma 2, we have, for ^ e K, i = 1, A'', 

N N 

(11) 

i=i j=i 

where Mi = DiJ2j'=i{'^Aj^ > 0. Now we are able to estimate Q[(/)~] for 
^ > miiij £*,i = 1, A^' as in (7) 



AT N 

,,3 U 



3=1 3=1 

Jm 

JR 



-Mie-^Mi+'') f ki{y)e^^^ydy 
Jr 

Jr 

= (e)+e-^^=«(9z^;^,(l-e 



(12) 



/ h{y)e^^^ydy) 
Jr ' 



For ^ > mini4*j e^'"'"^)"^"'^ is bounded above. Finally, from (12) and the fact 
that #(7ylc) < c, we conclude that there exists g > 0, which is independent of 
^, such that, for ^ > ^* 

Qi^(r )](? + c) > i^Xe"''^^ - 9i^;^.e-^^=«. (13) 
And since 0r(^) = for ^ < i = 1, A'' 



This completes the proof. □ 
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4.2 Proof of Theorem 1 (iii) with monotonicity of F 

Theorems that guarantee the existence of traveling wave solutions for cooper- 
ative systems have been established (e.g. [21,33]). In this section, it is assumed 
that F is non-decreasing on [0, j3] and from this assumption, we proceed to 
establish Theorem 1. 

As we state in Section 2, even for the case of monotone operators, the results 
and analysis in this manuscript are different from [21,33]. Here, we are able 
to characterize explicitly the asymptotic behavior of traveling waves through 
a careful analysis of eigenvalues and upper-lower solutions ( an analysis not 
provided in [21,33]). As we shall see, the analysis of the asymptotic behavior 
of traveling wave solutions for monotone operator enable us also to prove the 
existence of traveling wave solutions for non monotone operators. 

In order to complete the last step, we need to make use of the following 
Banach space 

Bp = {u = K) : u' G C(M), sup|u*(0|e''^ < oo,z = l,...,iV}, 
equipped with the weighted norm 

N 

\\u\l = Y,^up\u\0\e'^, 

where C(R) denotes the set of all continuous functions on R, and where p is 
a positive constant such that p < Ac- It follows that e Bp and € Bp. 
Finally, the following set is required (domain of the operator of interest): 

Ap = {u:ue Bp,r{^) < u{0 < ^^iO,^ e K} 

It is clear that Ap C Cp. By the standard procedure (see [26,13,28]), it can be 
shown that Qc[F] is a continuous map of the bounded set Ap into a compact 
set. 

Lemma 4 Assume {HI) — {H3) hold. Then Qc[F] : Ap ^ Ap is continuous 
with the weighted norm \\.\\p and relatively compact in Bp. 

Now we are in a position to prove Theorem 1 when F is monotone. Define 
the following iteration 

ui = {u\) = Qc[F{cj)+% u^+i = «) = Qc[F{un)],n> 1. (14) 

From Lemmas 1, 3, and the fact that F is non-decreasing, u„ is non-increasing 
on R, it follows that 

c(0 < <+i(0 < <(e) < 0+(e),e e r, « > i,« = i,...,iv. 

By Lemma 4 and monotonicity of {un), there is u G Ap such that limy^_^QQ [j'^^n — 
= 0. Lemma 4 implies that Q[u] = u. Furthermore, u is non-increasing. 
It is clear that limj^oo u'(^) = 0,i = 1,...,N. Assume that limj_,._oo Wi(^) = 
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fc*, i = 1, N k'^ > 0,i = 1, N because of u G Ap. Applying the dominated 
convergence theorem, we get ki = fi{k)- By (H2), k = (5. Finally, note that 

We immediately obtain that 

lim uXOe"^"^ = , i = 1, AT. (15) 

This completes the proof of Theorem 1 when F is monotone. 
4.3 Proof of Theorem 1 (iii) 

We proceed to characterize traveling wave solutions when the assumption that 
F is monotone is dropped. Our treatment is different even for the scalar case 
(A^ = 1). The key mathematical ideas used can be found in the literature albeit 
there are differences. Our use of the Schauder Fixed Point Theorem and the 
construction of the bounded set Vp are different from those found in [13,22]. 
As in Section 4, both Qc[F'^] and Qc[F~] are monotone. Note that F. F+.F^ 
have the same linearization at the origin. In view of the results in Section 4, 
there exists a non-increasing fixed point U- = € C^- of Qcff'"] such that 

Q4F-{u-)]=u- 

and lim{_)._oo = A"'* = l)---)-^) and limj^.oo w!_ (0 = 0,i = 1,...,N. 

Furthermore, limj^tx> w!_(C)e^'=^ = i'^^, * = 1, N. Let 

0+(O = (0?), 

where 

</.+ = min{^+, i^y-^^^}, C G K, i = 1, N. 

According to Lemma 1, </>+ is an upper solution of Qc Also note that if /3+ 
is replaced with </'+(^) is an upper solution of Qc[F~]. By the construction 
of it then follows that 

u_(0<^(O.Cgk 

Now let 

Vp = {u:u= K) G Bp,uUO < < ^iO,^ G (-oo,oo),i = 1,...,N}, 

where Bp is defined in Section 4.2. It is clear that Vp is a bounded nonempty 
closed convex subset in Bp. Furthermore, we have, for any u = (u*) e V 

w_ = Qc[F-{u-)] < Qc[F-{u)] < Qc[F{u)] < Qc[F+{u)] < Qc[F+{^)] < ^. 

Therefore, Qc[F] -.Vp^Vp. Note that the proof of Lemmas 4 does not need 
the monotonicity of F~ . In the same way as in Lemmas 4, we can show that 
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Qc[F~] : Dp — Bp is continuous and maps bounded sets into compact sets. 
Therefore, the Schauder Fixed Point Theorem shows that the operator Qc[F\ 
has a fixed point u in Dp, which is a traveling wave solution of (3) for c > c*. 

Since u^{C) < < ^ {— 00,00), i = 1, A'', it is easy to see that 

for i = 1, A'', limj_^oo = 0> limj-j-oo u'{^)e'^'=^ = v\^, 

/3r < liminf < limsupu'(0 < Pt 
and < < < A^,^ e (-00,00). □ 



4.4 Proof of Theorem 1 (iv) 

Proof The proof in this subsection follows the approach found in [4, 13]. We 
make use of the results in Theorem 1 (iii). Hence, for each m G N, we choose 

Cm > c* such that limm^oo Cm = c* . According to Theorem 1 (iii), for each 
Cm there is a traveling wave solution Um = {Um) of (3) such that 

Um = Q[F{Um)]{(. + Cm)- 

and 

lim u\0 = 0, p- < liminf <(0 < limsupw^^lO < Pt = 1. 

By the standard procedure (sec [26,13,28]), (wm) is cquicontinuous and uni- 
formly bounded on K. Hence, the Ascoli's theorem implies that there is a 
vector valued continuous function u = (u*) on R and subsequence {umu) of 
{um), such that 

lim Mm JO = 

K—>00 

uniformly in ^ on any compact interval of M. Further, the use of the dominated 
convergence theorem guarantees that we have 

u= Q[F{uM + c*) 

Because of the translation invariancc of Um, we always can assume that the 
first component ul^{0) equals to a sufficiently small positive number a > 
for all m. Since there is only a finite number of equilibria, we can choose a 
in such a way that it is not the first component of any nontrivial cquilibriTim. 
Consequently u is a nonconstant traveling solution of (3) for c = c* . □ 
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4.5 Proof of Theorem 1 (v) 

The proof of this subsection follows the approach in [13,22]. Suppose, by con- 
tradiction, that for some c G (0, c*), (3) has a traveling wave Un{x) = u{x — cn) 
such that u € with liminfa,_,._oo ^(a;) ^ and m(-|-oo) = 0. Thus u{x) can 
be larger than a positive vector with arbitrary length. It follows from Theorem 
l(ii) 

liminf inf Un{x) > /3~, for < c < c* 

n^oo \x\<nc 

Let c e (c, c*) and x = cn. Then 

lim u((c — c)n) — hm u„(cn) > liminf inf Un{x) > P~ . 

n—^oo n— >oo n—>oo \x\<nc 

However, lim„^oo u{{c — c)n) = u{oo) = 0, which is a contradiction. □ 



5 Minimum speeds and traveling waves for a competition model 

Hassell and Comins' model of the growth and spread of two population den- 
sities at time n and location x under an interference competition regime is 
used to highlight the applicability of the results in this manuscript. Wc makes 
use of the local analysis results of their model reported in [11]. The addition 
of the possibility of dispersal via the re-distribution kernel ki{x — y) leads to 
(2). If their two densities arc denoted by Xn{x) and Yn{x) then their model 
is given locally by a set of nonlinear coupled difference equations (1) with the 
addition of dispersal leading to (2). The following results highlight the contri- 
butions that the main theorem makes towards increasing our understanding 
of the role of dispersal, in the context of local competitive systems. Li [23] also 
investigated the minimum speed of (2). 

Model (2) can support four constant equilibria: the unpopulated state 
(0,0); the second-species monoculture state (0,r2); the first monculture state 
(n, 0); and (i-^r^^, t^^^)- The change of variables p = X,q = r2-Y allows 
to convert system (2) into the following coupled system of integro-difference 
equations 



Pn+i{x)= / ki{x-y)f{pn{y),qn{y))dy 

qn+i{x) = / k2{x -y)g{pn{y),qn{y))dy. 
Jm. 



(16) 



where 

5(P,9)=r-2-(r2-(z)e«-'^^f 
h{p) = pe~^ 
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It is clear that (2) and (16) arc not monotone systems. A straightforward 
calculation shows that (16) has fom- equilibria (0, 0), (0, r2), (ri,r2) and 

ri - C7ir2 ri-crir2, 
("i >o-2Ti ). 

1 — (Tl(T2 1 — (Ti<T2 

In fact, under the conditions of Theorem 3, we show in Appendix that there are 

no positive equilibrium of (16) between (0,0) and (ri,r2). Theorem 1 is used 
to guarantee the existence of a spreading speed and traveling wave solutions 
of the nonmonotone system (16) with its accompanying results on the speed 
of propagation.. We summarize the results obtained in the context of this 
example in Theorem 3. Its proof is outlined in the Appendix. 

Theorem 3 Let < r2 < 1 < ri, < cti < 1 < (72, CTicr2 < 1, and 
and 

Assume that k\,k2 satisfy (HI) and J^ki{s)e^''ds > J^k2{s)e^^ds for iJ > 0. 
Then the conclusions of Theorem 1 hold for (16). 

The biological interpretation of the conditions in Theorem 3 in the context 
of our application are straightforward. For an invasion to be successful, the 
overall dispersal of the invader (X) is relatively larger than the overall disper- 
sal of the out-competed resident (Y). Further competition favors the invader 
whenever cri is sufficiently small (invader less affected by competition) and 
a2 is sufficiently large (a relatively fragile resident, that is, more susceptible 
to interference competition). Under these conditions, there are traveling wave 
solutions of (16) "loosely" connecting its two equilibria (0,0) and (ri,r2). 
Equivalently, there are traveling wave solutions of (2) "loosely" connecting its 
two boundary states (0,r2) and (ri,0). Here the term "loosely" means the 
traveling waves may oscillate around the equilibria since they are not neces- 
sarily monotone. For specific fc;. the exact value of c* can be computed and 
compared to experimental data as it has been done by Kot, Lewis, others and 
their collaborators. 

6 Conclusions 

integro- difference systems arise naturally in the study of the dispersal of pop- 
ulations, including interacting populations, composed of organisms that repro- 
duce locally via discrete generations and compete for resources, before dispers- 
ing . The brunt of the mathematical research has focused on the the study of 
the existence of traveling wave solutions and characterizations of the spread- 
ing speed in the context of cooperative systems. In this paper, we characterize 
the spreading speed for a large class of nan cooperative systems, formulated 
in terms of integro-difference equations, via the convergence of initial data 
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to wave solutions. The spreading speed is characterized as the slowest speed 
of a family of non-constant traveling wave solutions. The results are applied 
to the non-cooperative competitive system proposed by Hassell and Comins 
(1976) [11]. We are in the process of applying these results to additional eco- 
logical and epidemiological systems where the local dynamics are naturally 
non-cooperative with the hope that increasing our understanding of the role 
of dispersal in communities where the local dynamics are richer, more realistic, 
than those previously supported by the mathematical theory. 
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Appendix 

Proof of Proposition 1 (4). 

If ^(yu) = ^ In ^(a*) achieves its minimum at a finite fj,, then c* = minj(i>o ^(m) > 
0. Now let c* lim^^oo ^ In A(/i). We recall that A(/i) is an eigenvalue of i?^ 
with a positive eigenvector. Thus there exists a positive constant 6 > and a 
positive integer i < N such that A(/i) > 6 J^ki{x)e'^^dx. Thus 



also see Weinberger [34]. Note that if'(O) = and therefore, c* > lim^^oo '^(a*) > 
0. □ 

Proof of Theorem 3. 

We verify that the conditions (H1-H3) hold for (16). From the assumptions 
of Theorem 3, (HI) holds for (16). We proceed to verify (H2) for (16) which, 
as we had noticed earlier, has four equilibria (0,0), (0,r2), (ri,r2) and 




Let tp'(/x) 
c* > lim^. 




■oo 




( 



n — air2 



n - air2 
1 — <Jia2 



)■ 



(17) 



1 — a\a2 
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If it is further assumed that ri > 1, r2 < 1, and ai < 1,<J2 > 1,<jiG2 < 1 then 

ri-(Tir2 ^i-o-i^2x^. . .1„^ 
,cr2- )>(ri,r2). (18) 

i — tTilT2 i — OiG2 

Thus (16) has no other positive equilibrium {p,q) between (0,0) and (ri,r2) 
with p > and q > Q. Observe that 1 is the maximum point of h{p), that 
is, /i(p) is not monotone on [0,ri]. Further simple calculations show that 
9p{v,q) = a2(r2-(7)e9-"^f > 0, for q e [0, r2], 5g(p, </) = (l-r2+(z)e9-'^^P > 0. 
In order to use Theorem 1, we define the upper monotone function 

Kp), o<p<i, 

/i(l) = e-S l<p. 
and corresponding monotone systems with /i+ 



h*{j,) 



qn+i{x)= I k2{x-y)g{pn{y),qn{y))dy. 
Jr 



(19) 



where f+{p,q) = h+{p)e''^-'^^''^+''^o. 

The origin, (0,0) is an equilibrium of (19) and g{p,q) = q has only two 
possible solutions q* = r2 and q* = (72P* ■ Thus for q* = r2, Equation (19) 
has two equilibria (0, r2), (e''^"-^, r2). The second equilibrium (e''^"^, r2) comes 
from the fact that p* > 1 and therefore h+{p*) = e~^. (If < < 1, 
then /i+(p*) = h{p*) and {p*,q*) = (ri,r2); however, ri > 1, which is a 
contradiction). In order for Equation (19) to have another positive equilibrium 
{p*,q*), when q* = G2P* ■, it must satisfy p* > 1 (otherwise, p* < 1 and 
{p*,q*) is (17) which means that (18) implies that p* > 1, a contradiction) 
and therefore 

q = a-2P . 

We will use the inequality, > a; + 1, :r e R to estimate e^. Thus p* = 
eri-(Tir2+ai(72P*-i > _ ^^^2 + <JiCF2P* and p* > ^jE§^ > ri, which implies 
that p* > e^^~^, from the first equation of (20) and (T2?'i > ?'2- Again since 
(J2r\ > r2, we also have q* = a2P* > r2 and thus Equation (19) has no positive 
equilibrium between (0,0) and (e'"^~^,r2). 

There is a G (Oi 1) such that h{to) = /i(e'"^~^) and define 



h-ip) 



Kp), o<p<to, 

h{to), to<p<e^'-\ 
and corresponding lower monotone system 

Pn+i{x)= / ki{x-y)f~{pn{y),qn{y))dy 

Jj (21) 

qn+i{x)= / k2{x-y)g{pn{y),qn{y))dy. 
Jr 
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where f-{p,q) = /i-(p)e''i-'^i''=+<"i«. Then 

< h-ip) < h{p) < h+{p) < h'iO)p,p e (0,e''i-^] 

h~{0) = h'^{0) = 0, h^{p),h(p) have the same derivative at 0. 

Since g{p, q) = q has only two possible solutions q = r2 and q = (J2P and 
/i(to) = ^oR *° = h{cf^~^) ~ e''i~^e^'^ ^ , wc can therefore calculate that 
(21) has throe equilibria (0,0),(0,r2) and (ti,r2) where ti = g'^ri-i-e^'i ^ 
Again (fi,r2) comes from the fact that ti > to and h{ti) = h{to)- (The same 
argument applied to (19) implies that ti < to is a contradiction. ) We will now 
show that 

ti<n. (22) 

Indeed, since Expression (22) is equivalent to 2ri — 1 — e^^~^ < Inn, we let 
l{x) = 2x — 1 — e^~^ — In a; and therefore = and 

l'(x) = 2 - e^-i - -<2-x--<0, for X > 1, 

X X 

and this verifies (22). Since < to < 1 < ri, the following inequality holds 

< e'^-'-'^'~' <to = e*°e'^^-ie-'="" < e^^^-'-'^'~' = h. (23) 

If {p*,q*) is a another positive equilibrium of (21) when q* ~ (^2P* , then it 
must satisfy p* > to (otherwise, p* < to and {p*,q*) is (17) and from (18) we 
have that p* > 1 > to, a contradiction) and therefore 

„ri-l -e^^^^ ri-mr2+cri<T2P* _ „* 

* (24) 

q = 0-2P • 

Since ri > air2 and p* > 0, we have p* > e^^~^~'^' ^ \ In view of the assump- 
tion, a2e'^'^~^~^ ^ > r2 and System (24), we have q* > r2- Again from (24), 
we have p* > e^^^~^~^ ^ = t\. Thus (21) has no other positive equilibrium 
between (0,0) and (ti,r2) and 

(0,0) « (ti,r2) < (ri,r2) < (e"-\r2), 

and this verifies (H2) for (16). 

We now proceed to verify (H3) for (16). The matrix in (4) for (16) is 

" -^"1^ >- \r2a2 /r k2{s)e^'ds (1 - rs) k2{s)e^''ds J ^^^^ 
Since e'^^~'^^^^ > 1 > 1 — r2, the principal eigenvalue for the matrix is 

A(At) = e"''-"'''' I ki{s)e^"'ds 
Jm. 
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and the corresponding positive eigenvector 

/,,(!) \ / e'-i-i'-2 J^fci(s)e^°ds-(l-r2)fBfc2(«)e'"ds \ 
7Ji,= i J2) 1 = ( r,a,f^k2{s)e^'ds 1 . (26) 

Because J^k2{s)e'^''ds < /j^ fci(s)e^''ds, canceling J^k2{s)e'^'^ds in vj^^ leads 
to 

,a)> ^""""-(i-^^) >l + ^""""-i >l (27) 

r2(T2 CT2 r2(T2 CT2 

It is clear now that (H3)(i) holds. We can proceed to verify (H3)(ii) for (19). 
Let 

(p, q) = (minje''^^^, i^j^^a}, niin{r2, a}), a > 0. 
Since e''~'^^'^ > 1 + 9 — <^2P, wo need to show that 

r2 - (r2 - (/)e«-"^P < raazp + (1 - r2)q + q{q - a2p) (28) 
< r202P+ (1 - r2)q 

Therefore, it is easy to see that we only need to verify that 

q < a2P (29) 

and 

< e-^i« (30) 

P 

For (29), we need to consider the two cases: p = e^^~^ and p = uj^^a. If 
p = e"~-^, then 

q<r2< aae'-i-i (31) 

which is true by the assumption (r2 < (T2e'^^~^~^'^^ '). If p = ujl^a, then 
q < a < a2i^^^a, which is true because of (27). 

In order to verify (30), first assume that p € (0,1), then h'^{p) = pe~^ 
and p = uj^^a since e'^^~^ > 1. Since e~'^^" < e~'^^'', it suffices to verify that 
g-i^f, < e"*^!", which is true because of (27) and (T21T1 < 1. For the case 
p > 1 we have h^(j>) = e"^. Again since 

it suffices to verify < < e~^^ , which holds because r2 < 1. 
It remains to veriiy (H3)(ii) for (21). Let 

ip, q) = (minjii, u'^^a}, min{r2, a}), a > 0. 
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For (29), wc need to consider the two cases: p = ti and p = f )i a. If p = ti, 
from the assumptions, we have 

q<r2< aae'^i-i-*'^'"' < ase^'-i-^-^"'"' = (Jaii = aap (33) 

and then (29) holds. If p = vjx^a, then q < a < a^i^j^^oe, which is true because 
of (27). 

We must verify (30) (with /i+ being replaced by h^) for (21). li < p < to, 
then h-{p) = pe'P and p = z/^^^a because of (23). Since e"''!" < e'"^", 
it suffices to verify that e""'' " < e"'^^", which is true because of (27) and 
a2<Ji < 1- For the case that p > to, then h~{p) = h{to)- From the definition of 
h~ and (32) we see that it suffices to verify that 

^LM < MM = e-*° < e--- (34) 
P - to - ^ ^ 

holds, which follows from Expression (23) and the assumption, 

e-'" < g-e^^-'-"'''' < e-''ir2_ (35) 

We observe here that the assumption air2 < e'^^~^~^''^ ^ can be relaxed as 
long as (34)(io > crir2) holds. 

To verify (H3)(iii), we note that h'^{p) ~ h^{p) — h{p) for p small, and 
conclude from Lemma 2, for sufficiently larger k, there is a small w » 0, 
if < {p,q) < u, f{p,q) > /(p,0) > (1 - i)e'-i-i'-^p and g{p,q) > (1 - 
i)r2CT2P+{l-i)il-r2)q. 

This concludes the proof of Theorem 3 since the conditions (H1-H3) have 
been verified. □ 
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